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This article gives a brief overview of some recent progress in the characterization and parametrization of density 
matrices of finite dimensional systems. We discuss in some detail the Bloch- vector and Jarlskog parametrizations 
and mention briefly the coset parametrization. As applications of the Bloch parametrization we discuss the 
trace invariants for the case of time dependent Hamiltonians and in some detail the dynamics of three-level 
systems. Furthermore, the Bloch vector of two-qubit systems as well as the use of the polarization operator 
basis is indicated. As the main application of the Jarlskog parametrization we construct density matrices for 
composite systems. In addition, some recent related articles are mentioned without further discussion. 

1. Introduction 
1.1. Motivation 

The state of a quantum system can be mathematically represented as a density matrix, a positive 
semidefinite hermitian operator of trace one. The density matrices of finite dimensional systems, 
which are the topic of this review, can be expressed as complex nxn matrices, constrained by the 
hermiticity, positivity, and trace conditions. Considering the fundamental role of density matrices 
for the description of physical systems, it is not surprising that many investigations have been 
devoted to the parametrizations of density matrices over the last years. This work is motivated by 
many reasons. For example, a suitably chosen parametrization may considerably simplify solving 
a specific physical problem, it may help to identify new properties of the system, or it can be used 
to study the properties of density operators itself. Furthermore, some parametrizations provide 
a straightforward way to generate positive matrices. When a general expression for a density 
matrix is needed, the positivity condition is typically the most difficult property to be sure of. 
Some parametrizations provide a simple way to overcome this problem. 

Maybe the best known density matrix parametrization is the Bloch vector parametrization. 
It was first used to describe the states of a two-level system, but was later extended to higher 
dimensions. An early review article discussing the properties of the Bloch vector was published 
by U. Fano in 1957 Q). In the standard approach using hermitian basis matrices, the Bloch 
vector of an n-level system is a real vector with — 1 components. In a two-level setting, the 
Bloch vector is a three component vector with length smaller or equal to one. It provides a way 
to identify the states of a two-level system with the points of a ball of radius one, the pure states 
corresponding to the surface of the ball. This mapping has found applications in hundreds, if not 
thousands, of articles. The two-level case will be mentioned in this review only very briefly as the 
main emphasis is on higher-dimensional extensions of the Bloch vector description. The structure 
of the set of Bloch vectors corresponding to positive operators becomes very complicated as soon 
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as the dimension increases from two to three. This is a consequence of the fact that when n > 2 
the maximal length of the Bloch vector producing a positive operator depends on the direction 
of the vector. The complicated nature of this set is visible already in the three- level case (d^. 
However, despite the problems associated with the structure of the set of Bloch vectors giving a 
physical state, the higher dimensional Bloch vector description has found applications in many 
fields, such as the dynamics of n-level systems. It has been shown that in this context the Bloch 
vector parametrization helps to identify constants of motions of some quantum optical systems 
Assuming that the Hamiltonian has a specific shape, the time evolution equation for the 
Bloch vector can be written in a block form. Vector components belonging to different blocks 
evolve independently of each other, and it turns out that the length of the vector inside each 
block is time independent (0, Q, [!)■ There exists also a class of constants of motion known as trace 
invariants. Unlike the constants of motion mentioned above, trace invariants are independent of 
the Hamiltonian of the system. The values of these invariants can be related to the components 
of the Bloch vector (Q). Another field where the Bloch vector parametrization has been used 
during the recent years is that of quantum entanglement. In particular, the entangled states of 
two qubits have been studied using this approach by various authors, see, for example, (Q, [l3"G3)- 
In the aforementioned applications the Bloch vector is defined in a basis consisting of hermitian 
operators. This is often a natural choice as it guarantees that the components of the Bloch vector 
are real. However, it may be preferable to use the polarization (or spherical tensor) operator basis 



nlyi . \2ui ] when the angular symmetries of states are important (l^il ). Another possible choice for 
the basis is the Weyl operator basis (f^ ). As in the case of the hermitian basis, defining the 
parameter set corresponding to physical states is a complicated problem also when the latter two 
bases are used. 

It is possible to express every n x n density matrix in terms of a diagonal matrix and an 
element of the unitary group U{n). This provides a way to parametrize the density matrices if a 
parametrization of U{n) is known. There exists many such parametrizations, but, depending on 
the on the problem, one of them may be preferable by providing a better insight or by leading 
to a reduction of the number of parameters. In general, a parametrization of U{n) has v? real 
parameters. Consequently, the size of the parameter set grows fast as n increases, encouraging to 
find ways to reduce the number of parameters. Two examples of parametrizations which allow to 



identify redundant parameters are the Jarlskog parametrization ^2B - \24 ) and the parametrization 
presented in (f^l ). These parametrizations have potential applications in quantum information 
theory, where many quantities and properties of quantum systems are obtained by optimizing 
the values of some functions over the set of all density matrices. This is the case, for example. 



when entanglement is quantified (|^a . \2'A \. Determining the values of entanglement measures 
leads often to numerically demanding calculations due to the high-dimensional parameter set 
over which the optimization has to be performed. Finding a parametrization which allows to 
construct a sufficiently large set of density operators while keeping the parameter set small 
would simplify his task considerably. The parametrizations discussed in this review are the coset 



^2a ) and Jarlskog parametrizations ^22424 ) ■ The former consists of parametrizing the density 
matrices in terms of certain cosets in the group U{n), while in the latter approach the elements 
of U{N) are given recursively, meaning that the elements of U{N) are expressed in terms of 
the elements of U{N — 1) and a unitary matrix containing the additional parameters needed to 
describe an element of U{N). With the help of the Jarlskog parametrization it is straightforward 
to generate matrices which are guaranteed to be positive. Furthermore, this parametrization can 
be easily extended to composite systems ([H). 

In addition to these approaches, there is a parametrization of (special) unitary matrices in 
terms of generalized Euler angles (3^, 31). It has been applied in [sF i to derive a volume formula 



for U{N) and related groups. This parametrization allows to eliminate redundant global phases 
in several cases. Yet another way to parametrize the elements of U(N) is in terms of quantum 



Householder reflections ^32 ) . This approach has been shown to be useful in quantum computation 



and quantum state manipulation ([5a45j). It seems, however, that these two parametrizations of 
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the unitary group have not been used to parametrize density matrices. Consequently they will 
not be discussed in more detail in this article. Another related factorization of unitary matrices 
in terms of orthogonal matrices is described in (fg^ ). A substantial generalization of these ideas 
can be found in (jga ). 

1.2. Basic definitions 

Recall that in quantum physics states are represented by density matrices on the complex Hilbert 
space Ti of the system. The Hilbert space Ti is always separable and in many cases of modern 
applications actually finite dimensional. 

A density matrix p on a Hilbert space Ti is by definition a linear operator from T-L into H such 
that 

p>0 i.e., {x,px) > for aU x e Ti (1) 

and 

oo 

Tr(p) = 5](e,-,pe,) = l, (2) 
j=i 

where (•, •) denotes the inner product of the Hilbert space Ti, {cj : j £ N} is any orthonormal 
basis of Ti and the normalization condition for the trace is written for the infinite dimensional 
case. In this case one has also to assume that p is of trace class. If Ti is of finite dimension n, 
then the sum in ([2]) extends from j = 1 to j = n. 

In this article we consider the case where the Hilbert space H is n-dimensional, n E N, 

U = U^ = C' . (3) 

We denote the basis of C" by {|1), |2), . . . , Now linear operators A 0x1% are represented by 
n X n matrices with complex entries, and in the following we do not distinguish between a matrix 
and the linear operator it represents. We denote hy Mn the space of all n x n matrices with 
complex entries. Then a density matrix p on %n is an element p £ Ain such that the positivity 
condition ([TJ and the normalization condition ([2|) for the trace hold. 

It is well known (and straightforward to prove) that in a complex Hilbert space a linear operator 
A which is positive in the sense of ([TJ (often also called positive semi-definite) is Hermitian, i.e., 
it satisfies A = A* , where A* is the adjoint matrix defined by 

{A*)ij = CLji for ah j = 1, . . . , n (4) 

Here aij are the coefficients of the matrix A and ciji denotes the complex conjugate of aji. 

The eigenvalues of a matrix ^ > (i.e., A satisfies ([T])) are calculated as follows (see, for 
instance. Theorem 25.1.1 of ©): ^he first eigenvalue is given by 

Ai = sup{(x. Ax) : X e Hn, \\x\\ = 1}, (5) 

and one proves easily from this definition that there is a vector ei G Tin, ||ei|| = 1 such that 
Aei = AiCi. The second eigenvalue then is ({ei}"*" denotes the orthogonal complement of {ei} 
in Tin) 

A2 = sup{{x,Ax) : X e {ei}"*" C Tin, \\x\\ = 1}, (6) 
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and there is 62 S {ei}"*", ||e2|| = 1 such that Ae2 = ^2^2- This procedure can be iterated and 
produces eigenvalues 

Ai,...,A„; Ai > A2 > ... > A„ > 0. (7) 

These n eigenvalues are not necessarily distinct; they occur in this list as many times as their 
multiplicity requires. For the eigenvalues of a matrix j4 > with Tr(j4) = 1 we know in addition 
to ^ that 

n 

Ai > and ^ Aj = 1. (8) 

j=i 

We know that the system of eigenvectors ei, . . . , of ^ is a complete orthogonal system in 'H.n- 
Therefore, the transition from the standard basis of to the basis of eigenvectors is effected 
by a unitary n x n matrix [/, and we arrive at the spectral representation of the matrix A: 

A = U*Dn{Xl,...,Xn)U, (9) 

with 



^rt(Al, . . . , A„) 



/Ai ••• 0\ 
A2 ••• 

'■■ '■■ '■■ A„_i 

V • • • XnJ 



(10) 



being the diagonal matrix of eigenvalues. 

Thus a density matrix can be characterized as follows: 

A matrix p e A4l\ {Ai^ denotes the space of all Hermitian elements of A^„) with coefficients pij e C, 
i, J = 1, . . . , 71 is a density matrix if, and only if, 

(a) p is positive in tlie sense of ([T|), i.e. all eigenvalues Xj are nonnegative, 

A,=A,(p)>0 for j = l,...,n; (11) 

(b) 

n 

Tr(p)-5]p,, = 1. (12) 

Let us introduce the set P„ of all n-dimensional density matrices, 

Vn = {peM'^: p satisfies ([II]) and ([l2])}. (13) 

In this article we study the general form of the elements of P^- The main difficulty is the positivity 
constraint p > 0. There are a number of criteria a matrix has to satisfy for positivity. Since they 
are typically given in terms of inequalities they do not provide much information about the 
concrete form of a density matrix. 

In this article we review recent studies on the structure and general form of a density matrix. 
Ideally one would like to have a parametrization of density matrices in the following sense: A 
parametrization of a density matrix p G "Dn means the following: 



(a) Specification of a parameter set Qn C where m depends on n, i.e., m = m{n); 
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(b) Specification of a one-to-one and onto map Fn : Qn — > T^n- 

Clearly, the case n = 1 is trivial: Di = {1}. Therefore in the following we assume n > 2. 

Unfortunately, this type of a parametrization of density matrices is not yet available for general 
n > 3. Only the case n = 2 is fully understood. We discuss here the partial results of the three 
main approaches to this problem: 

• the Bloch-vector parametrization; 

• the coset parametrization; 

• the Jarlskog parametrization. 



1.3. Overview 

Although the case n = 2 has been studied a lot in the literature, we start the following section on 
the Bloch-vector parametrization with a review of this case. This serves as a preparation for the 
discussion on the general case n > 3. This discussion is based mostly on (0; Q) ^^'^ it provides 
a reference point for the representations presented in sections [3] and [H where it turns out that 
all parametrizations agree for n = 2, but not for n > 2. We show that, despite the fact that the 
determination of the parameter set Q„ is a difficult task, the Bloch-vector parametrization can 
be used to identify various constants of motion or to study the properties of two-qubit states. 

Obviously, the spectral representation ([9]) can be used to obtain a parametrization of density 
matrices as soon as a parametrization of unitary matrices is known. Several parametrizations of 
unitary matrices can be found in the literature. We mention those which are of interest in connec- 
tion with our problem of parametrizing densit y ni atrices. One of them is the coset parametrization 
(see, for instance, ([H)). It has been used in (fgi ) as a starting point for an analysis of the space 
of all density matrices for n = 4. This analysis gives a description of the geometry of this state 
space in terms of flag manifolds. However, no parametrization of individual density matrices is 
offered. This idea has later been used in i2a ) to propose a parametrization of density matrices 
in terms of certain cosets in the group U{n). Section [3] contains a brief summary and discussion 
about the extent to which this parametrization can be considered a parametrization in the sense 
given above. 

The next section introduces the Jarlskog parametrization (see This parametrization is 

closely related to the coset parametrization and again starts from the spectral representation 
([9]) but uses the parametrization of U{n) presented by Jarlskog (see ([H, 0)) instead of cosets. 
The Jarlskog parametrization is recursive and thus allows a recursive parametrization of density 
matrices. We apply the Jarlskog parametrization to composite systems and obtain several quite 
interesting cases of a parametrization of density matrices for composite systems ([H). Finally, 
section [5] contains the concluding remarks. 



2. Bloch-vector parametrization 

In this section we first describe the Bloch vector of a two-level system, and then extend the 
discussion to n-dimensional systems. We show what is required from the Bloch vector of a three- 
level system for it to desribe a physical state. We illustrate the use of the Bloch vector in 
identifying constants of motion of dynamical systems, and then show how the Bloch vector 
approach can be used to study the separability of a two-qubit system. We also briefly describe an 
alternative basis for the Bloch vector parametrization, known as the polarization operator basis. 
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2.1. The case n — 2: The Bloch/Poincare sphere 
It is not difficult to see that the Pauh matrices 

(yx = (ri= (^J , ay = (72 = (^^ , a, = a-i=(^^ , (14) 

together with the identity matrix I2 form a basis of the real vector space Ml^oi Hermitian 2x2 

h 
2 



matrices. Hence every A ^ M2 can be represented as 



A = ao/2 + AifJi + A2cr2 + AsfJa, 

where oq and Aj, j = 1, 2, 3 are real numbers. Since the trace of the Pauli matrices vanishes, the 
coefficient oq is determined by Tr(^) = 2aQ. Hence, such a matrix belongs to T)2 if and only if it 
has a representation 

A ^^\^\ 1 / 1 + A3 Ai - iA2\ ;^ 

^=2^2+2LV. = 2(^A, + iA2 I-A3J' 

in which the coefficients \j are chosen such that all the eigenvalues of A are non-negative. 
The eigenvalues of A can be calculated by finding the roots of the characteristic polynomial 
det(x/2 — A). Using the abbreviation A^- = ^Aj we find for this polynomial (A = (Ai, A2, A3) G 

and \\? = T.U^]) 

det(x/2 -A) = {x- - iVp = - X + ^ - (16) 

The roots of this polynomial are 

xi = ^(l + |A|) and X2 = ^(1-|A|). (17) 

While the root xi is always > ^, the root X2 is non-negative if and only if |A| < 1. This gives our 
first result: The parameter set for this case is 

g2 = {A G : |A| < 1} = B(r3), (18) 

where i?(]R^) denotes the closed unit ball in with center at 0. We define a map on Q2 with 
values in P2 by the right hand side of (jlSp . i.e., 



I { 1 -I- A3 Ai - iA2 



By construction, this map is onto ^2- A simple argument shows that F2 is also one-to-one, i.e., 
i*2(A) = -F2(A') implies A = A'. We conclude that that {Q2,F2) is a parametrization of D2 with 
m = 3 = — 1 for n = 2. 

The parameter set Q2 = B{^) in this parametrization of P2 is called the Block or Poincare 
hall. In order to relate it to the standard representation and interpretation we just have to 
introduce spherical polar coordinates in B{ 



Ai = r sin 6* cos (/>, A2 = r sin 6* sin (/>, A3 = r cos 61, < 6* < vr, < </) < 27r, < r < 1. (20) 
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Using the spherical coordinates we can introduce a new parameter set 

^2 = {{e, (j),r) eR^ : 0<e<TT, < < 27r, < r < 1} (21) 
and on it the map F2 with values in T>2, defined by 

The boundary of the Bloch/Poincare ball is called the Block/ Poincare sphere. In terms of the 
parameter sets Q2 and Q2 this sphere is given by |A| = 1 and r = 1, respectively. Introducing the 
angle •& = 9/2 and using some elementary relations for the trigonometric functions sin and cos, 
the pure state density matrix corresponding to the point {9,(p, 1) on the Bloch/Poincare sphere 
has the form 

m^, <t^. 1) = ( "'^r^) = (23) 

where we used the abbreviations c = cost? and s = sin-!?, and 

1^) = (c,se'<^)^ = cosi?|0) +e''^sini?|l). (24) 

2.2. The general case n > 3 

As the last subsection shows, the Bloch vector parametrization works quite well in the case of 
2- level systems (n = 2). It has been generalized to the case of n-level systems with n > 3. 
The starting point of this generalization is the observation that the Pauli-matrices in (jlSp can 
be considered as the generators of the special unitary group SU{2). Therefore the following 
representation of density matrices p £ Dn has been suggested: 

1 1 

where Aj, j = 1, . . . , — 1, are the (orthogonal) generators of the special unitary group SU{n) 
and the Aj, j = 1, . . . — 1, are real numbers. The generators A-,- are n x n matrices with 
complex coefficients satisfying 

A* = Xj , TrXj = 0, Tr(A, Aj) = 26ij, i,j = l,...,n^-l (26) 

and the commutation and anti-commutation relations 

[Aj, Aj] = 2i ^ fijk^ki = —^ijln + 2 ^ Qijk^k- (27) 

fc=l fc=l 

Here the fijk and gijk are the structure constants of the Lie algebra su{n). It is also known that 
the generators Aj , i = 1 , . . . , — 1 together with the unit matrix /„ form an orthogonal basis 
of M.^ with respect to the Hilbert-Schmidt inner product {A,B) = Tic{A*B). This means that 
the inner product of two density matrices p and p, with Bloch vectors A and A, reads simply 
(p, p) = A • A 

The vector A is in this review called Bloch vector. In order to make a difference between the 
Bloch vector of a two-level system and that of an n-level system, n > 2, the term generalized 
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Block vector is sometimes used in the literature to describe the latter case. Yet another possible 
name for the Bloch vector is coherence vector. 

One possible explicit construction for the generators Xj can be given in terms of the generalized 
Gell-Mann matrices by defining the basis as in the following equation. 

\j){k\+\k){j\, l<j<k<n, 
-i\j){k\+i\k){jl l<j<k<n, 

+ + \ , 1<1<7^-1. (28) 



/(/ + ! 



If n = 2 these correspond to the Pauli matrices (|14p and if n = 3 they are the Gell-Mann 
matrices, see ()^9|) below. 

According to the properties (p6|) of the generators Xj every matrix of the form (j25p has trace 
1, Tip = 1. Thus a matrix of the form (|25p is a density matrix if and only if all its eigenvalues 
are non-negative. As in the case n = 2 these eigenvalues are the roots of the characteristic 
polynomial det(x/„, — p). As a polynomial of degree n the characteristic polynomial has a unique 
representation of the form 

n 

det(x/„ -p) = ^{-lyaj x"-^ ao = 1, (29) 

j=0 

where the coefficients aj are uniquely determined by the generators Xj and the parameters A = 
(Ai, . . . , An2-i) G M"^"^. In order to emphasize the dependence on A we write aj = aj{X). 
If xi, . . . ,x„ denote the roots of the characteristic polynomial det(x/„ — p), then (j29p shows 

n n 

^{-lya.x--^ = ll{x-x,) (30) 

j=0 j=l 

and by evaluating the product and comparing coefficients, the basic relation between the coeffi- 
cients aj and the roots Xj follows: 



a 



1 



Xi^Xi^ ■ ■ ■ Xi- Vieta's formula. (31) 



l<ii<i2<---ij 



It is known (an elementary proof is given in the appendix of Q)) that from this formula the 
important characterization of the non- negativity of the eigenvalues follows: 

Xj > 0, j = 1, ■ ■ ■ ,n, Oj > 0, j = 1, . . . ,n. (32) 

Accordingly we define the parameter set as Qn 

Q„ = {A G M"'-i : aj{X) > 0, j = 1, . . . , n}, (33) 

and on Qn the map Fn : Qn — ^ T^n by 



71^-1 
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The inverse of this map is given by 

Pn^{p) = (-^1 = Tr(pAi),... ,A„2_i = Tr(pA„,2„i)), p £ P„. 



(35) 



Thus, {QniFn) is a parametrization for n-level density matrices p € It is called the Bloch- 
vector parametrization. 

Now we discuss the difficulties with this parametrization for concrete applications. Although 
the parameter set Qn is uniquely specified, it is not easy to decide when a given point A G M" 
actually belongs to Qn ■ The origin of this difficulty is that firstly it is quite a complicated matter to 
explicitly calculate the coefficients aj for larger values of j and secondly it is not easy to determine 
explicitly the boundary of Qn, i.e., the set dQn = {A € M"^"^ : aj{X) = 0,j = 1, . . . — 1} 
since this amounts to solving polynomial equations in A of degree n, in — 1 variables. 

UsingLie algebra techniques, the polynomials aj{X), j = 1,2,3,4, have been calculated explic- 
itly in U) (see also (Q)). They read as follows: 



l!ai = 1, 
2!a2 

3!a3 

4!a4 



n-l 1 



n 



n^-l 



2n 



i,j,k=l 



(n-l)(n-2)(n-3) 3(n - 2)(n - 3) , ^ ,2 , 3(n 



ri.2-1 



re' 



4re 



n^-l 



n 



i,j,k=l 



(36) 



Already for re = 3 the parameter set has not been determined explicitly though 03 (A) is 
known. If in (f36|) we insert the explicit values of the structure constants gijk for the chosen set 
of generators of su{3), 03 is given by (recall that re^ — 1 is 8 for re = 3) 

3!a3(A) = ^ {s - 18|A|2 + 27A3(A1 + A^ - A^ - A?) - gVsxI 

+9\/3A8[2(A? + A^ + A^) - (A^ + A^ + A^ + A?)] + 54(AiA4A6 + A1A5A7 + A2A5A6 - A2A4A7)} 

(37) 

Thus it is obvious that the set 03 (A) > and with it the parameter set Q-^ has not been 
determined explicitly. However in (jj) some 2-dimensional sections have been calculated and 
represented graphicall}0. They show the very complicated nature of these sets. 

The Bloch vector parametrization has found applications in many different fields: It was used 
extensively by Harriman to study the properties of density operators (d, 043)- In particular, 
he discussed the structure of the set of physical states Qn (U)- This work was later extended, 
and partly reproduced, by other authors (14> yi- been useful in studying the dynamics 

and constants of motions of re-level systems 0); Markovian dynamics of decaying two- 

level systems ( ^ -■^). entropy production (0), characterization of the reachable sets for open 



systems driven by unitary control (1561 ). and unitary orbits in the set of density operators (jjil ) 



^We remark that equation (31) of (0) contains typing errors. The correct expression is given by I l37t above. 
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Furthermore, it has been used in studying the properties of two- and n- qubit systems ([J, 
Irh , the geometry of the states of a finite-dimensional systems (fjl ) , the classification of density 
matrices into different types (j53 . 15-^ ). detecting the entanglement properties of bipartite quantum 
states (fjl ). characterizing the structure of the state space of a two-qubit system (sdt i and the 
steady states of open quantum systems (jj'/l ). 
In the following we briefly discuss some of these applications. 



2.3. Block parametrization and trace invariants 

As we show next, n-level systems can be associated with constants of motion which are indepen- 
dent of the form of the Hamiltonian and depend only on the initial state of the system. These 
constants of motion are obtained by tracing over p", n = 1, 2, . . ., and are therefore called trace 
invariants. These invariants can be expressed using the Bloch vector and their values are also 
related to the coefficients Oj of the polynomial (j30p (0, Q) . The following discussion is based on 

Before introducing the trace invariants, we derive the differential equation that determines the 
time-evolution of the Bloch vector. The time-evolution of an n-level system can be expressed in 
terms of its density matrix p which satisfies the Liouville equation 



■ dpjt) 
dt 



[Hit), Pit)], 



(38) 



where the Hamiltonian H is in general time dependent and we have set h = 1. It has been shown 
in (fj^ ) (see also (5£, 6i)) that under certain conditions for the Hamiltoniai:H this equation can 
be solved in the form 



Pit) = UtpiO)U: (39) 

where Ut , t £ M, is a, family of unitary operators (which in general is neither a group nor a 
semi-group) and where p(0) is the initially given density matrix at time t = 0. This form of the 
time evolution guarantees that one stays in the space of density matrices. We write p as in ([25|) 

1 1 

P=-^n + -^XjXj. (40) 



Here and in what follows the time-dependence of the states and Hamiltonians is not always 
explicitly indicated. The Hamiltonian H can be expressed in a similar way 

n^ — l 

^=^In + lY.hAj, (41) 



where 

ho = TiH, hj = TiiXjH), j = l,...,r? -I. (42) 



Using ()40p . ()4ip and the Liouville equation, the following equation for the time-evolution of the 



'^Hii) is bounded and self-adjoint and the map t t— > Hit) is strongly continuous 
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Bloch vector can be obtained 

Yl fijkhjXk- (43) 
j,k=i 

The structure constant fij^ is antisymmetric, which guarantees that the length of A is time 
independent. In addition to this constant of motion, also other conserved quantities can be easily 
identified. We denote the rank of p by r. Then the density matrix can be associated with r 
constants of motion, as can be seen by writing the density matrix at instant t as 

p{t)=V{t)DV*{t), D = Dn{pi,p2,---,Pn). (44) 

where D is a diagonal matrix of the eigenvalues of p and y is a map F : M — )• U{n). The 
eigenvalues are time independent and V gives the time-evolution of the state. Now 

Tr(/(t)) = T,{D^) =p\ + pI + ... + pI (45) 

defines a constant of motion for each integer fc > 0. The number of non-zero eigenvalues is equal 
to the rank of p, so there are r independent constants of motion, given hyTi{p^), k = l,2,...,r. 
Due to the way they are defined, these constants are called trace invariants. Instead of using 
(I44p . it is possible to prove the existence of the trace invariants using (j43p and the properties of 
the generators {Xj} of SU{n). The values of the trace invariants can be expressed in terms of 
the Bloch vector A. For k = 1, 2, 3, 4 they are 

Trp = 1 

Mp') = z + lm' 



n 2' 

Tr(/>') = A + £:IAr + 7A-(A0A) 



3\ ^ ^ I \ 1 2 ■ ^ 

~^^2^^-^ ■ 4^ 
Tr(/) = + ^lAI^ + . (A A) + i^lAl^ + i|A Ap, (46) 



where we have defined 



{aQb)k = gijkaibj. (47) 



Explicit formulas for the trace invariants for n = 1, . . . , 9 given in terms of the Bloch vector can 
be found in the Appendix B of (Q)- Conditions for a vector A to describe a pure state can be 
obtained from these equations by setting Tt{p^) = 1. In this way one also obtains the boundary 
dQn of Qn, which consists of pure states. 

Comparing (j36p and (j46p we see that the coefficients can be given in terms of the trace 
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invariants. The expressions for the first four coefficients and the general expression are 



V.ai 


= Trp = 1 




2\a2 


= 1 - Tr(p2) 




Sla-s 


= 1 - 3Tr(p2) + 2Tr(/j3) 




4!a4 


= 1 - 6Tr(p2) + 8Tr(p3) - 


6Tr(/) + 3Tr(p2)2 


klak 




fc-i \ 
i=i / 



Clearly also the coefficients Oj are time independent. With the help of (|46p these can be expressed 
using the Bloch vector. 



2.4. Bloch parametrization and the dynamics of three-level systems 

We now describe how the Bloch vector parametrization can be used to study the dynamics of 
three-level systems. This topic was studied long before the structure of the set Q„,, that is, the 
structure of the set of Bloch vectors corresponding to physical states, was examined. We remark 
that in these studies it was not necessary to know the exact shape of the set Q„. This is a 
consequence of the unitary time-evolution, which guarantees that the Bloch vector X{t) belongs 
to the set Qn at all times t > if it is an element of Q„ at t = 0. One has only to make sure 
that the initial state determined by A(0) is a positive operator. 

Here and in what follows the initial time is chosen to be t = 0. The following discussion is based 
on • It was found in that the Bloch vector representation of the density matrices 

allows to identify, in addition to the trace invariants whose values do not depend on the shape of 
the Hamiltonian, also additional conserved quantities which exist only when the Hamiltonian has 
a specific form. We illustrate this by studying a three-level system, so the relevant Lie algebra is 
that of SU{3). One common realization for the basis elements {Ai, . . . , Ag} of this Lie algebra is 
given by the Gell-Mann matrices: 



Ai 





1 



1 



As = ^ 1 



0\ 




00 1 


1 


10^ 




/O 


' 







\0 i 






(49) 



The physical system we are interested in consists of a three-level atom interacting with two 
lasers. The energies of the three atomic levels are wi, W2 and and the energy differences needed 
now are denoted by uJi2 = oJi — UJ2 and u}23 = UJ2 — cos- The frequencies of the two lasers are 
vi and 1/2 and A12 = z^i — UJ12, A23 = 1^2 — ^23 are the detunings. We assume the case of exact 
two-photon resonance, which means that 



A12 = -A23 = A. 



(50) 
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13 



H 




(51) 



Here Q12 = • Si and = ^23 " ^2 are the Rabi frequencies, dj/^ is the atomic dipole moment 
between levels j and k and £1 and £2 are the vector amplitudes of electric fields of the lasers. We 
assume that the Rabi frequencies fii2, ^^23 have a specific form given by 



Qi2(t) = 2aOo(t) 
n23{t) = 2bno{t), 



(52) 
(53) 



where a and b are non-negative real numbers. This means that Q12 and $723 have the same time 
dependence but possibly different amplitudes. In the present case equation (|43p gives 



(54) 



where 



V 



1 








-A 


-b^Q 








\ 











-b^Q 





aOg 




















aOo 


A 








A 





















bn^ 

























— ailo 


-A 
























— 2aJ7o 













V 














-\/350o 





/ 



(55) 



and T denotes transpose. In order to simplify (154^ we define a new basis {A'^, . . . , Ag} through 
the equation 



(56) 



Here A = (Ai, . . . , Ag) and A = (A'j^, . . . , Ag) are vectors formed from the generators of SU{3) and 
the basis change is given by the time independent orthogonal matrix 



B 



+ 62 



fa 





b 


























a 





-b 












ab 












-1 






2a2+62 


V3b^ 





2Va^+b^ 







b 





—a 


























-b 





— a 











-a^+b^ 


















ab 


VSab 


Va-'+b^ 
\/3ab 


\/a^+b' 
V3b^ 


Va'^+b'^ 
-2a'^+b^ 




2Va^+b-^ 


2Va^+b^ 



(57) 



In the new basis the time-evolution can be solved from 



(58) 
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where V' is a block-diagonal matrix given by 

V' = BVB'^ 

and 




(59) 







-A 0^ 




^3' = 


A 


2e 








-2e ) 





/O -e \ 
e A 
-A -e 

\0 e J 



vi = 



(60) 



with 



(61) 



This result means that the time-evolution of the system can be analyzed in terms of three separate 
vectors: 

A3 = (A;,A'2,A'3), A4 = (Al,A'5,A'6,A'7), Ai = (A'8), (62) 
Because of the antisymmetry of the y -matrices, the lengths of these vectors are conserved 

IA4I = const., |A]^| = const. (63) 



I A3 1 = const.. 



These are not independent quantities as they are tied together by the normalization condition 
+ lAiP = lA'P =const. Therefore, in the presence of two monochromatic lasers and 



IA3P + IA 



ii4l 

under two-photon resonance, there are two conserved quantities in addition to the one arising 
from the conservation of |A'|. Due to unitary time-evolution A'(t) = A3(t) © A4(t) © Ai(t) assumes 
only values which produce a positive density matrix p{t). Using (|57p the conserved quantities 
(|63p can be expressed in terms of the components of the original Bloch vector A. For example, 
the time-independence of |Ai| means that 



const X (a^ -I- 6^)^, 



(64) 



where the value of the constant is determined by the initial values of A2, A7, and Ag. In ([3, 0) it 
has been shown that similar approach allows two identify two constants of motion in a three- level 
system under the assumption that the Rabi frequencies are sinusoidally modulated, with a phase 
difference 7r/2 between them. In this case the matrix determining the time-evolution is not block 
diagonal but enables nevertheless the time-evolution to be solved. 

This approach has been extended to systems with more than three levels (Q) ■ If the Hamiltonian 
has a certain form, obtained by defining a counterpart of the Hamiltonian (I5ip in the n level 
case, n conserved quantities can be identified. 



2.5. The Bloch vector of two-qubit system 

We now study the Bloch vector description of a system consisting of two qubits, denoted by 
A and B. The Hilbert space "H = TIa ®'Hb = (g) is four dimensional and the density 
matrices can be expressed using the the generalized 4x4 Gell-Mann matrices (|28p . However, 
instead of using this basis, it is often advantageous to define the basis in terms of the basis 
elements of the Lie algebra of SU{2) © SU{2). This approach has been used in many articles. 
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see, for example, ([J, u6l4-?a). The latter choice for the basis allows to examine the entanglement 
of two-qubit systems in a more natural way than the generalized Gell-Mann matrix basis. We 
choose the basis {Ai, A2, . . . , A15} of the Lie algebra of SU{2) (g) SU{2) as 

Ai = ^cJi®/2, i = 1,2,3 (65) 
v2 

Ai = -^h <S) o-i-3 i = 4, 5, 6 (66) 

Ai = -^ai ai-Q i = 7, 8, 9 (67) 
v2 

Ai = ^cj2 cJi_9 i = 10, 11, 12 (68) 
v2 

Ai = ^CT3 (g) f7i_i2 i = 13, 14, 15. (69) 
v2 

This forms an orthogonal basis with respect to trace with the normalization given by 

Tr(A^Aj) = 2 5ij. (70) 
The totally symmetric structure constants are 

5147 = 5158 = 5169 = 524(10) = 525(11) = 526(12) = 534(13) = 535(14) = 536(15) 

= -57(11)(15) = 57(12)(14) = 58(10)(15) = -58(12)(13) = -59(10)(14) = 59(11)(13) = (71) 

These are needed in the calculation of the coefficients of the characteristic polynomial. An arbi- 
trary two-qubit state can be expressed as 

1 1 

= + - J] Ai Ai (72) 

i=l 



^ 3 3 3 3 

—7= (^-^4 + ^ Ai (Ji (g) /2 + X] -^«+3 h®CFi + ^^ Aj+3i+3 di ® fJj^ , 
i=l i=l i=l 7=1 



(73) 



where A is assumed to be such that p is positive. The reduced single-particle density matrices 
determined by this state are 

1 1 , ^ 1 1 - ^ 

PA = TrB(p) = -/2 + ^ J]] AiO-i, PS =TrA(p) = -/2 + ^X^Ai+3ai (74) 

^ i=l ^ 1=1 

Here Tr^ denotes trace over subsystem A, Tib is defined similarly. The components of the Bloch 
vectors of the reduced states are uniquely obtained from the components of the Bloch vector of 
p. The inverse does not hold; any information regarding the Ay, As, . . . , A15 components of the 
Bloch vector is missing from the reduces states pj^ and p^- The knowledge of the Bloch vectors 
of PA and pb does no allow to construct a unique two-qubit Bloch vector which gives rise to pA 
and Pb- 

As an example of the use of the parametrization we consider the Werner state for two qubits 

I — X 

Pw{x) = — — h + xS, (75) 
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where x is a real parameter and S is 



s = i 

2 



/O 0\ 
1-10 
0-110 

\0 0/ 



(76) 



We first determine the range of the parameter x which ensures positivity of pw(x), making sure 
that pw{x) indeed is a density operator. After this we study the separability of the Werner state 
W{x). We may write Vl^(x) as 

1 X 

W{x) = -I4 - - (CTI (g) (71 + (72 (T2 + 0-3 (g) (73) 

/l-x \ 
_ 1 1 + x -2x 
~ 4 -2x 1 + x 

\0 l-x J 



The only non-zero components of the Bloch vector are Ay = An = A15 = - 
ward calculation gives the coefficients of the characteristic polynomial ([36 

V.ai = 1 

2[a2 = l{l-x') 



3!a3 
4!a4 



;(1 - Sx^ + 2x^ 



32 



(1 - 6x^ + 8x^ 



3x' 



(77) 

-x/^/2. A straightfor- 

(78) 
(79) 

(80) 

(81) 



These are all non- negative when —1/3 < x < 1, which is therefore the range of x corresponding 
to physical states. Usually, however, the range of x is taken to be [0, 1]. We now examine the 
separability of W{x). The density matrix p of a composite system (g is called separable 
when it can be written as a probabilistic mixture of tensor product states 



^PiPA ® Pb, Pi > 0, ^Pi = 1- 



(82) 



If a state is not separable it is entangled. Detecting separability is in general a very complicated 
problem, but in some special cases it can be done easily. A simple way to test the separability 
of a two-qubit system is to use the Peres-Horodecki, or positive partial transposition, criterion 
(6£, [g!l ). It provides a necessary and sufficient condition for separability for 2x2 and 2x3 
dimensional systems. According to this criterion, a state p is separable if the operator obtained 
by transposing the density operator of the subsystem A or B is a. positive operator. When the 
state of the second subsystem is transposed, (i73l) becomes 



[-T^h + '^Kcri^ 12 + "^ Xi+3 ^2 ® 0-f + ^ ^ Aj+3i+3 CJi ® a 
i=l i=l i=l j=l 



(83) 
(84) 



i=l 



i=l j=l 
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where 



-A,, A 



pt 



-A, 



.All 



All , A 



pt 

14 



-Ai4 and A?* = Aj for the rest of the compo- 



nents. When partial transpose is appHed to the Werner state we get 



Wptix) 



1 X 



(Ti - (72 (8) <T2 + 0-3 Cgl CT3) 



/I 



V 







2x \ 





1+x 
1+x 
-2x 1 - 2; 



(85) 



The calculation of the coefficients of the characteristic polynomial of Wpt{x) shows that a 
and are the same than for W{x) and 



pt 



3!af = -(1-3x2-2x3) 



4!, 



32 



(1 - _ g^3 _ 3^4^^ 



(86) 
(87) 



These differ from 03 and 04 by the sign of the x^ -term. These coefficients are non-negative 
when —1 < X < 1/3. Therefore the Werner state W{x) is separable when —1/3 < x < 1/3 and 
entangled when 1/3 < x < 1. Because W{x) is obtained by partially transposing Wpt{x), the 
state Wpt(x) is separable when —1/3 < x < 1/3 and entangled when — 1 < x < —1/3. 



2.6. The polarization operator basis 

Our definition of the Bloch vector with respect to the generalized Gell-Mann matrices is not the 
only possible. Although changing the basis will change the Bloch vector representation of the 
states, it does not affect the structure of the set of density matrices in general. One possible way 
to choose the basis is to use the polarization operators, also known as spherical tensor operators. 
The concept of a polarization operator appears in the quantum mechanical th eory of angular 



momentum and is thoroughly discussed in the literature, see, for example, (jg^ . \6a ). When the 
angular symmetries of the system are important it is convenient to expand p using the polarization 
operators. Many examples of such systems can be found in [2! ). We will describe the polarization 
operators only briefly. In the following discussion we define the basis and characterize some of 
the properties of the set of Bloch vectors corresponding to physical states. This subsection is 
based on ([13, [13). 

The polarization operators pertaining to an n-level system are defined as 



Tlm = xl^^ECZ:,LM\k){l\, (88) 
k,i=i 



where the indices have the properties 

n — 1 



2 ' 

L = 0,1, ... ,2s, 

M = -L,-L + 1,...,L-1,L, 

mi = s, 7712 = s — 1, . . . , nin = —s . (89) 
The coefficients C^*"*" r are identified with the usual Clebsch-Gordan coefficients C;]"! • ^ of 
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the theory angular momentum. For L = AI = the polarization operator is proportional to the 
identity matrix 

Too = -^In. (90) 



All other polarization operators are traceless, but they are not Hermitian in general. Due to the 
symmetry properties of the Clebsch-Gordan coefficients they satisfy the orthogonality relation 

^"^{Ti^MTL^hh) = h^LjhhKh (91) 

and they have the property 

tIm = {-^)''Tl-m. (92) 

The former equation ensures that the polarization operators form an orthonormal basis. Any 
density matrix can be written using this basis as 

P=l^- + flY. ^lIiTlm = + AP° • T (93) 

L=l M=-L 

with the Bloch vector in the polarization operator basis given by A''° = 
(A^*" ^, A^Q, A^°, A2°2) ^2° 1' -^20' -^('n-i)(n-i))' '^^lerc the components are ordered and given by 
'^LM ~ 1'^(-^1m'°)' general, the components X^j^j are complex since the polarization operators 
are not Hermitian. The hermiticity of the density matrix, p = p* , forces the components of the 
Bloch vector to fulfill the condition 

AfM = (-l)''Ar_M (94) 

In particular, the components A^q are real. In order for p to describe a physical state it also has 
to be positive. As in the case of the Gell-Mann matrix basis, the positivity of p can be checked 
using (129p and (132^ . The coefficients aj of the characteristic polynomial (129p are related to the 
trace invariants through the equation 

ja, = ^(-l)— ia,_^Tr(p'"). (95) 

m=l 

General expression for the trace invariants can be calculated to be 

1 , fc(fe-l),,p„,2 , ^ /A;\Tr[(AP°-r) 



j^k—m 



m=3 ^ ^ 
n— 1 Li n—1 Lm 

Tr[(AP°.r)-]= E •••E E ^i:M.---^ZM„J^iTL.M.---TL^Mj. (96) 

Li = l Mi=-Li L„=1M,„=-L„ 



The traces can be calculated using the equations presented in (jga ). Combining the results we 
find that the coefficients of the characteristic polynomial are 



k=2 



^ + E 



"A:\T¥[(AP°-T)™] 



1 Z-^ \m) m 

m=2 



(97) 
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The first three coefficients are ao = oi = 1 and 



2a2 = 1 |AP°|^ 

n 



(98) 



Therefore a necessary condition for p to be a positive operator is that Bloch vectors he within a 
hypersphere of radius \X^°\ < \/ {n — l)/n. In the case of a two-level system (n = 2) the density 
matrix becomes 



p = -I2 + (a + i/3)rii - (a - i/3)Ti„i + jTqo, 



(99) 



where a,/3,7 are real and we have defined = a + if3, Aqq = 7 and used The set of 

physical states is now given by Q2 = {(a,/3,7) £ \2{a^ + (3^) + 'y^ < 1/2}. The surface of this 
set corresponds to pure states and is a prolate spheroid. 

As in the case of the Bloch vector given in the generalized Gell-Mann basis, for n > 3 the 
structure of the set of physical states becomes very complicated (f7^ ). Nevertheless, pure states 
are on the surface, mixed ones lie within the volume and the maximally mixed state corresponds 
to |AP°| = 0. 



3. The coset parametrization 



This section presents a short summary of the article n28i ). If D = Dn{Xi, . . . , A.„) is the diagonal 
matrix of eigenvalues in the spectral representation ([9]) of a density matrix p, we denote by D' 
the commutant of D in U{n), i.e., the set of all matrices U € U{n) such that UD = DU. Now, 
if ([9]) holds for some U G U{n) and V is some other unitary matrix such that VU~^ G D' then 
one has 



p = U*DU = V*DV, VU~^ £ D', 

and conversely. 

In the case of nondegenerate spectrum the commutant D' is easily determined: 



(100) 



D' = T" = C/(l)(8)---(g)f7(l), n factors, (101) 

and thus one can say 

p = n*Dn, n£U{n)/T''. (102) 

According to the book (fg3 ) elements U € U{n) can be factored in the following way: 

u = nA~i ■ ■ ■ ^2^1 (103) 

with $7i G T" and 

U{k) T'^-'^ 



Typical coset representatives Vt}, are of the form 

' SU{k)/U{k-l) 



k = 2, . . . ,n . 



(104) 



^n—k 



(105) 
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where O is the k x [n — k) zero matrix and its transpose while In-k denotes the n — k 
dimensional unit matrix. 

After a brief sketch of the general case the article (f^3 ) proceeds to discuss in some detail the 
cases n = 2 and n = 3 (and the Bures metric and the n = 2 state space). A more systematic 
approach and an explicit realization of the factorization (jl03p is given in the following section. 
Therefore no further details are presented here. 



4. The Jarlskog parametrization 

Like the coset parametrization also the Jarlskog parametrization starts from the spectral represen- 
tation ([9]) of a density matrix. Accordingly, we begin with a concrete version of the factorization 
of a unitary matrix ?7„ G SU{n) in terms of certain basic unitary matrices as in (1103^ following 
(0). Since this factorization is based on the use of canonical coordinates of second kind in the 
Lie group SU{n), a more accurate name for the parametrization would be the "canonical coor- 
dinate parametrization". Some preliminary investigations into this problem have been given in 
the article (j^ ). After having defined the Jarlskog parametrization of unitary matrices, we show 
how it can be used to parametrize density matrices. As an application of this parametrization 
we construct a general density matrix of a two-level system. We also show that the Jarlskog 
parametrization can be straightforwardly extended to composite systems and illustrate this by 
constructing some two-qubit states. 



4.1. Jarlskog 's recursive parametrization of unitary matrices 

Recall that the Lie group U (n) is connected but not simply connected and that every Un € U (n) 
has a determinant of absolute value 1. It follows that every Un G U{n) has a determinant 
det Un = and thus can be written as 



where SU{n) denotes the Lie group of unitary matrices of determinant 1. Obviously, the diagonal 
matrices D„(e*", !,...,!) and Z)„(Ai, . . . , A„) commute and hence in ([9]) we can restrict ourselves 
to unitary matrices in SU{n), i.e., /?„ G if) and only if, 



where the eigenvalues Xj satisfy (jT]) and ([8]). This reduces our problem to that of finding a 
parametrization of unitary matrices of determinant 1. 

In order to provide the necessary background for this parametrization of unitary matrices and 
the important recursion formula we follow Fujii (fg3 ) to explain the origin of the basic building 
blocks of this parametrization. He observed that the Jarlskog parametrization of unitary matrices 
is obtained by using canonical coordinates of the second kind for the Lie group SU (n) (see, for 
instance, (0)). 

Recall that SU{n) is the (simply) connected component of the unit element of U{n) and thus 
is the image of its Lie algebra su{n) under the exponential map. The Lie algebra su{n) consists 



C/„ = Z?„(e*",l,. 



,1)K, aGM, KGSUin) 



Pn = KDniXl, . . . , Xn)Un, Un G SU{n). 



(106) 
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of all skew-adjoint n x n matrices X. Such matrices have the form 





/ 

1 




Zl2 


Zl3 


Zl,n—1 


Zln \ 








-Zl2 


ia2 


Z23 ■ ■ 


• Z2,n-1 


Z2n 








-Zl3 


-Z23 


ia^ ■ ■ 


• 23,n.-l 


Z2n 


(107) 


X = 
















-Z\^n-\ 


-Z2,n- 


1 Zs^n-l • • 


• ian-1 


Zn—l,n 






I 


-Z\,n 


-Z2,n 


-Z3,n ■ ■ 


Zn—l,n 







where the a^'s are real numbers and the Zjk are complex numbers which can be chosen indepen- 
dently. Such a matrix has the natural decomposition 



X = Xi + X2 + • • • + Xj- + • • • + X„ 

where X\ is the diagonal matrix with the diagonal entries ia\ , . . . , ia^ 

X\ = Dn{iai, . . . ,ian) 



(108) 



(109) 



and where for j = 2, . . . , n the matrix Xj is the matrix which has in column j the column vector 



/ Zij \ 

Z2j 



:iio) 



and in row j the row vector 



{Zj\ = {-Zlj,-Z2j,...,-Zj-lj) 



;iii) 



as entries. All other entries of Xj are zero. 

The canonical coordinates of the second kind for SU (n) are given by 

su{n) B X = Xi + X2 + ■ ■ ■ + Xj + ■ ■ ■ + Xn ^ e^'e^' ■ ■ ■ e^' ■ ■ ■ e^" £ SU{n) 

The exponential e^^ is easily calculated: 

A,i = e^^ = A,i(ai,...,a„) = Z)„(e^°S...,e^<^"), G M. 

For j = 2, . . . , n — 1 we write 



(112) 



(113) 



X, 



"n-j 



\zj) 
-{Zj\ O 



(114) 



where Ofc denotes the kx k zero matrix and O indicates that the remaining entries of the matrix 
are zero. Naturally Xn = K^- Thus we get 



, j = 2,...,n-l, 



(115) 
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where denotes the k x k unit matrix. Using 



^1 



(\^^^^^^\ ^ - -(z-\z-)K- 



one can calculate the exponentials and finds for j = 2, . . . ,n, \zj\ = ^ {zj\zj), 

1 „9 . , , 1 



e 



= Ij + (1 - cos |— i^i^l + sin \zj\--Kj . (116) 



If one introduces the unit vector \zj) = jj-^\zj) G ^ one can rewrite this exponential as 

e^^=I,-{l-cos\z,\)(\'f^\ j)+sin|z,|(^|- Ij^) . (117) 
The right hand side of (jll7p is a unitary matrix Vnj which is often written as 



Vn,j = C'~' i"'^], c,=cos\z,\, s,=sin|z,|. (118) 



Since the adjoint (and the inverse) of the matrix e^^J is e~^^ , the adjoint and the inverse of the 
matrix Vnj is given by (|118p with Zj replaced by —Zj. 

The Jarlskog matrices Anj = e^' , j = 1, . . . ,n (see {2£, 0)) thus are of the form 



-^n,j — ^ Qj''' J ^ ) j — 2, . . . , n 1, ^n,n — Vn^n- (H9) 

A generic element G ^C/ (n) therefore has the factorization 

Un = e^-e^- • • • e^" = ^„,l^n,2 • • • ^n,n (120) 

or equivalently 

Un = e "e " ^ • • • e ^ = j4„^„j4n • • • A^^i- (121) 

Observe that by construction the Jarlskog matrix is defined in terms of the following set of 
parameters 

Oj > 0, Zj e s{a-^) 

where we changed notation: 9j for \zj\ and Zj for Zj, and where S{C^~^) denotes the unit sphere 
in We indicate this by writing Anj = Anj{6j, Zj). The special structure of the matrix Anj 

and the properties of the trigonometric functions imply that on the parameter set 

P^ j = {{ej,Zj) : < % < 7r/2, zj G S{a~^)} (122) 

the mapping {6j,Zj) — > Anj is injective, i.e., if {6j , zj), (6^, z^) G Pn.j and Anj{Oj, Zj) = 
Anji9'j,z'j) then Oj = O'j and zj = z'- 

According to (|107p or (|120p . (|12ip and (jl22p this factorization describes a generic element 
Un G SU{n) in terms of real parameters. 
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Notice that the matrices Vnj do not depend exphcitly on the dimension n and that we can 
write, for j = 2, . . . ,n — 1, 

I O 1 j' - I O In-j-lJ- ^ ^' 



It follows that 



^M---^n,.-i=(^"'^--/"-"-^V)- (124) 



The analysis presented above shows that Vn.i • • • Vn,n-i, up to a factor in T"~^, gives the generic 
factorization of (n — 1) x (n — 1) special unitary matrices and hence, from (|12ip . we get the 
important recursion relation 

TT -/^^"-I'^^/l _ fUn^l 0\ fin-l - {I - Cn)\Zn){Zn\ Sn\Zn)\ 

Next we address the question of injectivity of the Jarlskog parametrization of Un G SU{n) 
modulo elements in T", i.e., the parametrization (I120p without the factor A^^i G T". Then this 
parametrization uses the parameter set 

Pn = {{Oj, Zj) : {Oj, Zj) G Pn,j, j = 2,...,n} (126) 

where we have taken (jl20p and (jl22p into account. In order to prove injectivity of the map 

(6*2, Z2; Z3;...; 9n, Zn) > Un = ^n,2(6'2, 22)^n,3(^3, Z^) ■ ■ ■ An ,n{9m Zn 

where {62, Z2', ^3, z^; . . . ; On, Zn) G -Pn, we proceed by induction with respect to the order n. Since 
(9j,Zj) — 7- Anj is injective on Pnj, the map (jl27p is injective for n = 2. Now assume that for 
some n > 2 this map is injective for the orders k < n — 1. We now show that the map (I127p then 
is injective for the order k = n. 

Suppose that for Un, Un G SU{n) we have U'n = Un where U'n = Un{02, z'2, 0'^, z'^,; . . . ; O'n, z'^). 
In the recursion formula ()125p for C/„ we abbreviate In-i — (1 — Cn)\zn){zn\ with Bn-i and 
similarly for U^ and B'^_^. If we calculate the matrix product in the recursion relation, the 
identity {/^ = C/„ reads 



?7/,_iS^„l s'nU'^_^]\z[^)\ ^ fUn^lBn-1 SnUn-l\Zn) 
~^n{^n\ ) \ ~Sn{Zn\ Cn 



(128) 



where we used the abbreviations s'n = sin 9'n and = cos O'n- Thus = c„; since On, O'n G [0, 7r/2] 
we conclude O'n = On and therefore s'n = Sn- The identity —s'n{z'n\ = —Sn{zn\ now implies 
{z'n\ = {zn\ and \z'n) = \zn)- Next we use the identity s'nU',^_i\z'^) = SnUn-i\zn) to conclude 
U'^_i\z'^) = C/^_]^|2;„) = Un^i\zn)- Finally we use the identity U^^^B'^^i = Un-iBn~i to get 

U'n-l - (1 - Cn)U'n-i\Zn){Zn\ = Un-l - (1 - Cn)Un-l\Zn) {Zn\ 

and from the identities established above we find U'^_-i = Un-i - The induction hypothesis implies 
{O'j, z'j) = {Oj, Zj) for J = 2, . . . , n — 1, and we conclude. 
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4.2. Jarlskog parametrization 

Since diagonal matrices commute, (jl06p implies the following parametrization of a generic density 
matrix p„ G P„: 

Pn = U*Dn{Xl, \n)Un = ■ ■ ■ A\^Dn{\l, . . . , A„)A„,2 • • • (129) 

where the matrices A^^j and the matrix D„ are parametrized as described above. Recall that 
for each j = 2, . . . , n both the parameter set for the matrix An^j and the concrete form 
of this matrix have been given explicitly. The number of parameters is easily calculated: There 
are r? — n real parameters for the product An^2 ■ ■ ■ ^n,n and n — 1 parameters representing the 
eigenvalues Ai,...,An subject to the normalization condition X]j=i -^i ~ This gives ii^ — 1 
independent real parameters as in the Bloch vector parametrization. 
If we introduce the set 

n 

An = {A = (Ai,...,A„)GM": Ai > A2 > • • • > A, > 0, J]] A, = 1} (130) 

we can specify the parameter set Qn of this parametrization as 

Qn = {(A, Zj) : A G A„, £ [0, 7r/2], zj G 3(0-^), j = 2, . . . , n}. (131) 

Our recursion relation (|125p for unitary matrices leads to an interesting recursion relation for 
density matrices in different dimensions. Denote the product An^2 ■ ■ ■ ^n,n-i by Un-i- Then, using 
()125p . we can deduce from (jl29p the recursion formula 

_ (U*^,0\ (Dn-i{\i,...,\n-l) 0\ (Un-lO 

Pn ^n,n Q 1 ^ 1^ O A,-J 1^ O 1 

= K,n (^^Q ^ ^^^^^ 

where 

Pn-l = K.iDn^l{\l,...An^l)Un-l (134) 

is a positive matrix in dimension n — 1 with trace Trp„_i = 1 — A^. 

While in the case of the Bloch vector parametrization the question of uniqueness was quite easy 
to answer, it is fairly complicated in the case of the Jarlskog parametrization and only partial 
answers are known. Given A G A„, consider the diagonal matrix Dn{\) of eigenvalues and the 
commutant of this matrix in SU{n), 

Dn{\)' = {VCi SU{n) : VDn{\) = D.n{\)V}. (135) 

Clearly this commutant depends on A G A„. As mentioned in the previous section, in the case 
of a non-degenerate spectrum, i.e., if Ai > A2 > • • • > A„ > 0, then this commutant is easily 
determined 

n 

D„(A)' = {y = Z)„(ei"S...,e^"") : 5^0^=0}. (136) 




Next consider the case that for some 1 < k <n one has Ai > A2 > • • • Ajt > while Aj = for 
J = A; + 1,. ..,?!. This case can be reduced to the previous one by restricting the density matrix 
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to the subspace spanned by the first k eigenvectors. 

In general one can not exclude degeneracy in the spectrum of a density matrix p„,; then one 



has to consider eigenvalues Xj with multiplicity rrij, j = 1,2, . . . , k, k < n, Xi > X2 > 
mi + • • • + TTT-fc = n. For this kind of a spectrum the commutant (jl35p is 



X, 



Dn{Xy = {V 



V2 



VkJ 



: Vj £ U{mj), det y = JJ det Vj = 1}. 



(137) 



In order to establish an injective parametrization of density matrices in this case one has to 
determine the Jarlskog parametrization of elements C/„ € SU{n) modulo elements in Dn{Xy . 
This is not yet known. 



4.3. Simple examples 

As has been explained in section 2, the structure of the set of Bloch vectors is complicated when 
n > 2. In order to overcome this other parametrizations can be used. To get some inspiration 
on how to proceed we have a new look at the Bloch-vector parametrization for n = 2. In the 
parametrization (j22p we again take 9 = 2^0 and write the relations (|17p for the eigenvalues as 
|A| = xi — X2 and 1 = xi + X2- Then, with the abbreviations c = cosi? and s = sini9, the 
parametrization (|22p can be written as 



xic^ + X2s'^ {xi — X2)cse " 

{xi — X2)cse^'^ X2C^ + XiS 



,2 , . .2 h (xi,X2,i?,.^) GQ2 (138) 



where the parameter set Q2 has the form 

Q2 = {(xi,X2, !?,(/>) G : < xi,0 < X2, xi+X2 = l, < i9 < 7r/2, < (/>< 27r}. (139) 
In this form we can recognize the parametrization (|138p to be the following product of 5 matrices: 
^e-W2 \ (c -s\ (xi ^ \ ( c s\ (e"^/"^ , 



e'<^/V c ) \ Q X2) \-s c) \ e-'"t>''^ 

This is of the form p = U*D2{Xi, X2)U where the diagonal matrix of the eigenvalues .02(^1, A2) 
is sandwiched between a unitary matrix U and its adjoint U* when we set Xi = Aj. This agrees 
with the parametrization (|129p for n = 2 with 



U = A 



2,2 



1 - (1 - C2)|2;2)(^;2| 52^2) 

-S2{Z2\ C2 



when we set Z2 G ^(C) as Z2 = e''^ and 6*2 = G [0, 7r/2]. 

Knowing the explicit form of a density matrix p2 in two dimension we can easily calculate the 
explicit form of a density matrix on by using the recursion formula (I132p : Given any point 
{>^iA2,h]62,Z2;03,Z3) G Q3 we find 

P3 = ^3,3(^3, ^3)*^3,2(e2, ^2)*^3(Al, A2, A3)A3,2(^2, ^2)^3,3(^3, ^3) (Ml) 
= ^3.ft °>3,3, (142) 
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where p'2 is a positive matrix with trace Trp2 = Ai + A2. 

Clearly, this recursive construction can be continued to n = 4, 5, , 



4.4. Jarlskog parametrization of density matrices for composite systems 

Consider now a density operator for the composite system living in ® C". It is clear that 
we may parametrize it as a density operator living in C"™. However, by doing this we lose 
information about the particular tensor product structure of the total Hilbert space C^*". To 
control the division into subsystems of C"™' = ® C" let us consider p as an n x n matrix with 
m X m blocks, i.e.. 



Pn,m = ^ ® Pij, 



(143) 



with pij being m x m complex matrices. Our aim is to provide a suitable parametrization for 
positive block matrices. We proceed in analogy to the previous section. 

For A = (Ai, A2, . . . , Xnm) £ ^nm wc denote by Dnm{K) the diagonal matrix of eigenvalues. 
Then for any Un^m £ SU{nm) 



Pn,m — Un,m-^ni'n{^Un^-, 



(144) 



is a density matrix in C"™. Again we parametrize elements Un,m G SU{nm) by 



su{nm) 3 Xi+X2 + ---+Xn — ^ e^^e^' • • • e^" G SU{nm 



(145) 



where now Xi is an nx anti-hermitian block-diagonal matrix with m x m blocks, and for j 
2, . . . , n is an nx anti-hermitian matrix with m x m blocks defined as follows: 



J I -^m 



(146) 



where instead of the n — 1 column vectors Zj used in the previous section we choose now n — 1 
column block vectors 



(147) 



with Zij being m x m matrices; and similarly 
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The Jarlskog matrices An,m = are now of the following form: 

/Ui Om-- - Om\ 

^mU2 ■■■ Om 



A 



where the Uk are unitary m x m matrices. For j = 2, . . . ,n one finds 



27 



(148) 



/vl 



\ 



(149) 



where the m dimensional unit matrix appears n — j times and where Vn,m is a unitary j x j 
block matrix with m x m blocks, defined as follows: 



' n,m 



Here we use the following notation: \Zj) denotes the normalized block vector 

- 1 2 * 

I > \\Zj\\ = ZljZij -\ h Z*_ 1 J Zj-ij, 



(150) 



z, 



(151) 



and (S) Im — \Zj)[Ij-i CS" {Im ~ Cj)]{Zj\ stands for the {j — 1) x (j — 1) block matrix 



\k){l\®Zlfi,Zi, 

k,l=l 



31 



(152) 



with 



C j — cos S j — sm J , — II Zj 1 1 . 
Our parametrization (|144p of density matrices Pn,m now reads 

— 4"* ... /1 2* n fA, I ... I A )A^ • • • 4" 

where 

Dn{^l \ ■ ■ ■ |A„) = Al^ ,|^Dnm{K)^n,m 

is a positive block diagonal matrix with diagonal blocks 



(153) 



(154) 



(155) 



(156) 
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which are positive m x m matrices and which satisfy the normahzation condition 

Tr(Ai + ---A„) = 1. 



(157) 



4-4- 1- 2x2 systems 

Let us consider a 2 Cg) 2 system to illustrate our parametrization for the well-known 2-qubit 
states. Taking 



Ai = O2, A2 = ^{h - CTz), S = smal2, C = cosal2, U = az 
one obtains a family of rank-1 projectors 



(158) 



^ sin^ a sin a cos a\ 

sin a cos a 

P{a) = 



ysin a cos a cos^ a j 



(159) 



which corresponds to a pure state 

^pa = sin a 1 00) + cos a |11). 

Note that this state is separable if and only if S = or C = 0. For S = C = , one obtains 

a maximally entangled state. It shows that a nontrivial rotation by a does produce quantum 
entanglement. 

As a second example in this class let us take S = C = I2/ V^, U = ax and 



Ai 



I fl-p 

1 - W ' 



Ao 



I fl-p 
l + 3p 



(160) 



with —1/3 < p < 1 to guarantee positivity of the matrices Aj. One obtains the partially trans- 
posed Werner state (I85p 



Wpti-p) 



fl+p 2p \ 

l-p 

l-p 

\ 2p l+pj 



(161) 



As has been shown earlier, this state is separable if and only if p < 1/3. The point p = 1/3 is 
not distinguished by our parametrization. 

Next we consider S = sin a I2 and C = cos a I2 and obtain a more general two-parameter 
family 



/(p, a) = -h 0I2+ pP{a) 



(162) 



of 2 (g) 2 states. This family is separable if and only if 



1 

- 1-F2sin(2a)' 
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This family of states can be generalized further as follows: Instead of (jl60p we take 



Ai 



P2 



A, 



P3 



(163) 



0p47' " \0 pi 

where pt >0 and pi + P2 + P3 + Pi = ^- Furthermore, the S and the C matrices are chosen as 

sin a \ ^ I cos a 



S 



sin /3 y ' 



c 



cos /? 



a,^ G [0,7r/2]. 



With U = ax the following family of states results (sq, = sin a, Cq, = cos a and similarly for 
and Ci3, p = {pi,P2,P3,Pa) as above): 



a, 13) 











/ Picl+P2s'i 


(P3-P4)sQ;Ca P3S^p+PiC^p 

\{Pl- P2)siiCp 



(Pl -P2)Sf}Cp\ 




Pl4+P24 / 



(164) 



By construction or by a direct calculation we see that p{p] a, /3) > and Trp(p; a, /5) = 1 for any 
choice of the parameters with the restrictions given above. 

Note that the above family belongs to the class of 2 (g) 2 circulant states considered in (fg^ . 
Note also that for a = /3 = 7r/4 the family of states (jl64p reduces to the family of Bell diagonal 
states 



(Pl +P2 



pip) = 2 



Pl- P2\ 



P3+Pi P3-P4 

P3 - P4 ^3 + ^4 
\Pl -P2 pi+ P2J 



(165) 



In (|g<j ) the separability of these states has been investigated for various values of the parameters. 

4-4-^- 2 (S> m systems 

For n = 2 and m > 3 our formula (jl54p reads 



/1 2* /Ai\/i2 /i2 2 / ZCZ* ZS 

P2,m = -42,m-^2(Al|A2)A2,m> ^2,m = ^2,»i 



(166) 



with Z = U £ U{m) and C = C0SH2 and S = sinH2. -D2(Ai|A2) is given by (llSSp . (jl56p . and 
(|157p . Matrix multiplication gives 



P2,i 



U* 



CU*AiUC + SA2S SA2C -CU*AiUS^ 
CA2S - SU*AiUC CA2C + Su*AiuS J 



U 



(167) 



By choosing particular values for the parameters appearing in (|167p we find some examples which 
have been considered in the literature. 

For S" = O or C = O one obtains a class of block diagonal matrices 



A 



1 

'm 



A2 



or 



A. 



2 

K Ai 



:i68) 
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These matrices represent separable 2(g)m states and thus show that quantum entanglement arises 
only for nontrivial H2 such that C 7^ O and S ^ O. 

Next consider the case Ai = A2 = A and [A, U] = 0. We get the following class of 2(g) m states: 



U* Orn] (A B\(U Om 
Im) \B* A) \Om Im 



(169) 



with 

A = CAC + SAS, B = SAC = CAS. 
If in addition U*AU = ^ is assumed one gets the matrices 

iiUB)' ™) • 

These are block Toeplitz positive matrices and it is well-known that they are separable (fgj ). 
Thus (I170p defines a huge family of bipartite separable states. 

Similarly, the block Hankel positive matrices of (fgj ) can be reconstructed. To this end we 
assume that the matrices U, Ai , A2 and H2 satisfy 

[U*AiU,E2]=0, [A2,H2]=0. 

This produces the following class of 20 m states: 



U* Om\ (Ai B'\ ( U O 

Im) \B' A2 [Orn /, 



(171) 



with Ai = CAiC + SAiS and B' = SC {A2 - U*AiU). Under the additional assumption UB' = 
B'U we arrive at matrices of the form 

with X = UB' . These are block Hankel positive matrices and hence separable (fg^). 



5. Conclusion 

In this review we have discussed different parametrizations of nxn density matrices. We have com- 
pared three different parametrizations, namely the Bloch vector, coset, and Jarlskog parametriza- 
tions. Of these the Bloch vector parametrization is the oldest and most widely used, while the 
two others are relatively recent discoveries and consequently not so well known. 

The Bloch vector parametrization has found more applications in physical problems than any 
other density matrix parametrization. This is due to the inherent simplicity of the parametriza- 
tion. The basis matrices used in this representation are hermitian, which guarantees that the 
components of the Bloch vector are real. Furthermore, these components can be obtained straight- 
forwardly as expectation values of hermitian operators. Therefore, given a density operator, the 
calculation of the components of the corresponding Bloch vector can be done straightforwardly. 
On the other hand, there is a serious disadvantage in the Bloch vector parametrization: In all 
the cases of n > 3 it is practically impossible to determine the parameter set corresponding 
to physical states explicitly (see Section 2.2). Only for a two-level system the parameter set is 
easily determined and is given by unit ball in M'^, the so called Bloch ball. The difficulty derives 
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from the requirement that a density matrix has to be a positive operator. As yet, there is no 
simple way to determine which Bloch vectors lead to positive matrices. Despite the complexity of 
the parameter set, the Bloch vector parametrization has many applications. It has been known 
for a long time that it gives a bijective mapping between the states of a two- level system and 
the points of the Bloch ball, providing an elegant way to visualize the states and dynamics of 
two-level systems. Another advantage is that the hermiticity of the basis enables to write the 
density matrices and Hamiltonian in the same basis. This makes possible to derive easily the 
differential equation giving the time-evolution of the Bloch vector. Under some conditions for 
the Hamiltonian of the system, this approach allows to identify various constants of motion. 
The Bloch vector has turned out to be useful also when the entanglement and separability of 
two-qubit states is examined. 

The second parametrization discussed in detail in this article is the Jarlskog parametrization 
which is based on a suitable parametrization of (special) unitary matrices. Its main advantages 
are (see sections 4.2 and 5.1): 

• Its parameter set is given explicitly; 

• it is recursive; 

• it extends naturally to composite systems. 

However, in the form presented here, it contains redundancy in the case of a degenerate spectrum. 
We have indicated how to eliminate this and get injectivity for this parametrization, too. 
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